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LIST  OP  SYMBOLS 


-a 


X  coordinate  of  doublet  or  doublet  line 


a^  polynomial  coefficients  of  source  distribution 
b  radius  of  a  bulb  (approximated  as  a  sphere) 
d  length  of  source  line 

P„j  P,,  P^  Proude  numbers  with  respect  to  f^  d,  and  L  re¬ 
spectively 

f  depth  of  a  point  source  below  the  free  surface 
f^/^  depths  of  end  points  of  source  line 

g  acceleration  of  gravity 

G  (k  )jl  Integrals  defined  by  Equations  [l6]  and  [l4]  re- 
°  spectively 

K  jK  Modified  Bessel  Functions  of  the  2nd  kind 
0  1 

k^  =  g/V^ 

L  distance  between  two  vertical  source  lines  fore  and  aft 


m 


m  total  strength  of  a  point  source  or  a  source  line 
= 

r  radius  of  half  body 

R  total  wave  resistance 


R  ,  R  ,  R, 

0  s  b 

Rj  4.)  K  ■  J  R « 

int  £i  is 


R^  (a,b,  c),  R^(a,b,c,  )  j 
Re  =  real  part  of 


,  Parts  of  Wave  resistance  defined  in  Sec- 
!  tions  7  and  lO. 
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R1  [ak^jk^jn]  an  Integral  defined  by  Equation  [12] 

V  Uniform  free  stream  velocity  at  oo 

yi,Y,z  Rectangular  right  handed  coordinate  system  with  origin  at 
free  surface,  z  positive  upward,  and  x  In  the  direction  of 
the  uniform  flow  velocity  V 

<t)  total  velocity  potential 

4)^ , (t>2, 4)^0  potentials  defined  In  Section  2 

p  strength  of  doublet 

linear  coefficients  of  strength  of  doublet  In 

V  =  vyv 

fictitious  frictional  force 
p  mass  density  of  water 

m  =  (x  +  a)  cos  0  +  y  sin  a 

^  wave  height 

4  , L  wave  height  due  to  sources  and  doublets  respectively. 

S  U 
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ANALYSIS  OF  THE  BULBOUS  BOW  ON  SIMPLE  SHIPS 

1.  INTRODUCTION 

The  effect  of  bulbous  bows  on  the  wave  resistance  of 
ships  was  first  Investigated  by  the  systematic  experiments  of 
D.  W.  Taylor  (I9II  and  1943) j  E.  M,  Bragg  (I930),  and  E.  P. 
Eggart  (1935) .  It  had  been  generally  understood  that  the  de¬ 
crease  of  resistance  due  to  a  bulbous  bow  Is  a  wave-making 
phenomenon.  J.  G.  Thews  In  I93O-I932  conducted  experiments  to 
determine  the  conditions  for  the  cancellation  of  the  bow  wave 
by  the  bulb.  (See  H.  E.  Saunders  1957)-  In  1928j  Havelock 
calculated  the  wave  form  due  to  a  doublet  Immersed  In  a  uni¬ 
form  stream^  and  found  that  a  wave  trough  was  formed  just  aft 
of  the  doublet .  Since  a  deeply  Immersed  sphere  Is  equivalent 
to  a  doublet,  W.  C.  S.  Wlgley  (1936)  Investigated  this  effect 
with  a  mathematical  model  of  a  Mlchell's  ship  plus  a  doublet. 

By  using  Havelock's  resistance  formula  (1934b),  he  demonstrated 
the  fact  that  the  doublet  wave  cancels  the  bow  wave  and  thus 
lessens  the  wave  resistance.  His  analytic  work  was  also  sup¬ 
plemented  and  confirmed  by  his  model  experiments.  G.  Welnblum 
(1935)  dealt  with  this  problem  by  expressing  the  form  of  a 
ship  with  a  bulbous  bow  In  terms  of  a  polynomial  according  to 
Mlchell's  thin  ship  approximation. 

Recently  Inul,  Takahel  and  Kumano  (i960)  observed  wave 
profiles  and  found  that  the  wave  due  to  a  ship  with  a  submerged 
sphere  faired  Into  the  bow  was  exactly  the  superposltlo  .  of  the 
wave  due  to  the  corresponding  point  doublet  and  that  due  to  the 
hull.  They  explained  the  effect  of  the  bulb  on  the  wave  resis¬ 
tance  of  a  ship  by  using  the  Idea  of  Havelock  's  elementary 
surface  wave  (1934a). 
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The  mathematical  expression  of  the  wave  resistance  of 
a  conventional  ship  is  extremely  complicated.  Except  for  the 
special  cases  of  very  low  or  very  high  Froude  numbers  it 
is  necessary  to  integrate  a  highly  oscillating  function  numeri¬ 
cally.  And  it  is  extremely  difficult  to  investigate  the  effect 
of  varying  parameters  in  the  wave  resistance  equation.  There¬ 
fore,  only  simple  models  of  ships  are  considered  in  the  present 
paper.  The  advantages  of  this  procedure  are:  1)  We  can  analyze 
not  only  the  wave  resistance  but  also  the  effect  of  the  size 
and  the  location  of  the  bulb.  2)  We  need  not  be  concerned 
about  the  effect  of  linearizing  the  boundary  condition  on  the 
ship  surface.  3)  We  can  investigate  the  fundamental  relation¬ 
ship  between  the  source  and  doublet  under  the  free  surface. 

At  first,  the  simplest  system,  a  point  doublet  and  a 
point  source,  is  considered.  The  so-called  Interference  term 
of  the  wave  resistance  is  calculated  by  a  series  expansion. 

The  optimum  distance  betvifeen  two  singularities,  and  the  optimum 
size  and  depth  of  the  bulb  are  obtained.  By  this  procedure  we 
can  show  that  a  remarkable  reduction  in  the  wave  resistance  can 
be  realized  by  the  use  of  bulb. 

Second,  a  system  consisting  of  point  doublet  and  a  finite 
source  line  is  considered.  This  is  treated  in  similar  manner  to 
the  first  case . 

Third,  a  system  consisting  of  a  vertical  doublet  line 
under  the  free  surface  and  a  vertical  source  line  from  the  free 
surface  is  considered.  The  strength  of  the  doublet  line  is  con¬ 
sidered  to  vary  linearly.  However,  the  optimum  distribution  of 
the  strength  of  the  doublet  line  is  found  to  be  almost  uniform. 
The  Influence  of  a  stern  is  also  considered  using  the  method  of 
stationary  phase.  In  this  case,  the  calculation  is  performed  on 
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the  I.B.M.  1620  Digital  Computer. 

To  obtain  the  wave  resistance,  Lagally's  theorem  is  used 
In  the  first  case  while  Havelock's  formula  Is  used  in  the  second 
and  third  cases.  When  Lagally's  theorem  Is  used  we  can  see  very 
clearly  how  a  thrust  force  is  applied  at  the  point  of  the  doublet, 
when  It  Is  located  at  the  proper  position  In  front  of  the  source. 

The  actual  shapes  of  the  low  resistance  systems  are  found 
by  computing  the  actual  stream  lines.  The  bulb  Is  found  to  be 
very  large  so  that  the  bulb  Itself  may  be  considered  to  be  a  bow 
rather  than  an  appendage . 

This  method  of  analysis  is  shown  to  be  applicable  to  the 
general  case  In  which  the  waterline  of  a  ship  is  expressed  as  a 
polynomial . 

Case  1 

2.  FORMULATION  OF  THE  BOUNDARY  VALUE  PROBLEM 

The  origin  0  is  placed  on  the  mean  free  surface .  The  x 
direction  Is  the  same  as  that  of  the  velocity  of  the  uniform 
stream.  The  z  axis  Is  directed  upward  from  the  surface.  The 
coordinate  axes  0-xyz  shown  In  Figure  1  form  a  right  handed  sys¬ 
tem.  Symbols  are  shown  In  the  figures  or  In  the  symbol  table  If 
not  mentioned  In  the  text.  The  water  Is,  as  usual,  considered  to 
oe  incompressible,  homogeneous,  and  invlscid.  The  motion  Is 
considered  to  be  steady. 

A  system  consisting  of  a  doublet  at  point  (-a,0,-f),  a 
source  at  point  (0,0, -f)  and  a  sink  at  point  (L,0,-f)  In  a  uni¬ 
form  stream  will  represent  approximately  the  combination  of  a 
bulb  and  a  Ranklne  ovoid.  The  wave  resistance  due  to  this  system 
can  be  obtained  by  using  Lagally's  theorem.  To  Investigate  mainly 
the  relation  between  the  bulb  and  bow,  we  only  need  consider  the 
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system  of  the  doublet  at  point  (-a,0,-f)  and  the  source  at  point 
(0,0j-f).  We  denote  and  <t>^  as  the  perturbation  potential 

due  to  the  point  doublet  and  the  point  source  respectively.  The 
symbols  m  and  -p  denote  the  strength  of  the  source  and  the  doublet 
respectively . 

Denoting  R(ajbjc)  as  the  x  component  of  the  force  at  the 
point  (a^bjc)  we  obtain  by  Lagally's  theorem  (Lagally  1922  or 
see  Milne  Thomson  I956.)  . 


R^(-a,0,-f) 


=  -  4  TT  p  p 


Bx 


(-a,0,-f) 


[1] 


R  (0.0, -f) 

2 


4  TT  p  m 


B<t>  d<l> 
_ ^  20 

Bx  Bx 


(0,0. -f) 


[2] 


The  subscript  0  to  0^  means  that  we  exclude  the  effect  of  the 
point  Itself. 

Now  we  have  only  find  the  perturbed  velocities  at  two 
points  (-a,0,-f)  and  (0,0,-f).  The  perturbation  potential  0 
must  satisfy  Laplace's  equation. 


0  =  0 


[3] 


and  the  boundary  conditions  on  the  free  surface  and  at  infinity 
as  well  as  near  the  singularities .  The  linearized  kinematic 
boundary  condition  on  the  free  surface  [see  Lamb  1945]  is  given 
by 


0  =  -V  t  on  z=0 

z  X 

where  ^  is  the  free  surface  elevation.  If  the  idea  of  fict¬ 
itious  frictional  force  [see  Lamb  1945  or  Lunde  1952]  is  used, 


[4] 
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the  pressure  condition  on  the  free  surface  is 

4)V-gC  +  v<l>  =  0onz  =  0  [5] 

X  ^ 

V  being  the  coefficient  of  the  fictitious  frictional  force, 

1 

Combining  above  the  two  conditions  we  obtain 

(t*  +k  <t)  +  v<t>  =0  [6] 

XX  o  z  X 

on  z  =  0 

where  =  g/V^ ,  v  =  v^/V, 


When  the  distance  from  the  singularities  goes  to  infinity, 

V  <t>  =  0  [7] 

In  addition,  we  note  the  empirical  fact  that  the  disturbance 
in  front  of  ship  decreases  very  rapidly  when  -x  becomes  large. 


3.  FORCE  AT  DOUBLET  POINT 


The  solution  of  our  problem  [3]  with  boundary  conditions 
[6]  and  [7]  can  be  found  in  many  papers  (e.g,  Lunde  I952)  .  Us¬ 
ing  the  Integral  representation 


m 


V (x+a)^  +  y^  + 


=  —  / 
0 


2v 


d0 


I 

0 


-k(  z-f  -leu) 


dk 


where  cu  =  (x+a)  cos  0  +  y  sin  0,  we  have  for  the  potential, 
due  to  only  the  Images  of  the  doublet,  satisfying  the  conditions 
[3],  [6]  and  [7] 
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M.(x+a) 

"r  ^  ^  1  3/2 


^  TT  00  k(lCU-  Z-f  )  ,, 

„  o  r  r  1  k  sec  6  e  '  '  dkdS 

-Re  —  (J-  /  I - 

-IT  0  k-k  sec^  e  -  Iv  sec  9 
o 


where  Re  means  "real  part  of".  Hence 


(x+a)^  +  +  (z-f)‘ 


3  (x+a)‘ 


(x+a)"^  +  y^  +  (z-f)‘ 


k  TT  «>  ,  2  k(la)-  z-f  I)  ,, 

„  0  .  r  k  e  '  I  -^'dkde 

+  Re  —  /  /  - 

-V  0  k-k  sec  0  -  Iv  sec  0 

o 


The  X  component  of  the  perturbed  velocity  due  to  the  source, 
which  satisfies  Equations  [3],  [6]  and  [7]  Is 


1  X 

x,0,-f  ~  2  +  1 

""  x2+  4  f^i 


k  m  TT  00  .  ,  _ 

_  0  r  r  Ik  sec  0  e 
-Re  —  f  f  - — 


k(lx  cos  0  -  2f) 


dkd0 


-TT  0  k-k  sec^  0  -  Iv  sec  0 
o 


HYDRONAUTICS j  incorporated 


-7- 


Hence  from  Equation  [1] 


R  (-a^Oj-f)  =  -  4  TT  p  p. 


k  T  “  ,  ,  3  „  -2  k  f 

_  Or  r  ikcosSe  dkd0 

Re  q  —  /  /  - — - 

-TT  0  k-k  sec^  0  -  Iv  sec  0 
0 


+  m 


+ 


3a‘ 


a' 


a^  +  4f- 


3/2 


a^+4f^ 


5/2 


+  Re 


k  m 
o 

TT 


V 


/ 

-TT 


00 


/ 

0 


,  2  -k(la  cos  0  +  2f ) 
k  e  ^ _ Mkd0 

k-k  sec^  0  -  Iv  sec  0 
0 


The  integral  v;lth  respect  to  k  can  be  evaluated  using  the  contour 
Integral  as  in  Equation  [11].  After  the  integration^  v  is  set 
equal  to  zero.  Then  the  first  integral  can  be  represented  by 
means  of  modified  Bessel's  functions  (see  Lunde  1952)  as  in 
Equation  [15]  i.e. 


TT 


Re  / 

-IT 


00 


/ 

0 


—Pic  f* 

i  k^  cos  0  e _ dkd0 

k-k  sec“  0  -  Iv  sec  0 
o 


TT 


-k  f 
0 
e 


X 


K  (k  f)  +  1  + 


2k  f 
o 


5 


4 .  NONDIMENS lONAL  FORM 

We  may  use  the  relation  between  q  and  the  radius  of  the 
sphere  b  in  a  uniform  stream  at  infinite  depth 


[9] 
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Slmilarly 


r 


2 


where  r  Is  the  maximum  radius  of  a  half  body  produced  by  a 
source  In  a  uniform  stream. 

Then  In  the  nondlmenslonal  form  as 


R 

1 


F  =  =  Froude  number,  etc  . 


k 

0 


r 


r 


f 


and  dropping  bars  for  convenience  we  have 


R  (-a, 0,-1) 


2  1  3a^ 

+  i~7r 


1  1 

la^  +  4) 

a^  +  4 

+ 


2b‘ 


-k 


®  2 
F^r^ 
1 


K  (k  )  + 
o  o 


1  + 


2  k 


h  O'O’' 


-  Re 


2b" 


TT  F^ 


where  K  and  K  represent  the  zero  order  and  first  order 
o  1 

modified  Bessel  functions  of  the  2nd  kind  respectively,  and 


[10] 
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R 

0 


TT 

=  Re  / 


-TT 


«>  ,2  -k(la  cos  6  +  2f) 

k  e  '  dkd0 


0  k-k  sec^  e  -  Iv  sec  6 
o 


[11] 


5 .  EVALUATION  OP  INTEGRAL  R 

o 


Taking  a  contour  as  shown  In  Figure  2  for  the  Integra¬ 
tion  of  (l)  with  respect  to  k  and,  noting  that  v>0  for  the  resi¬ 
due  at  the  singularity  k  =  k^  sec^  9  +  v  sec  9  and  later  let¬ 


ting  v-k).  We  obtain 
tt/s 


R  =  2Tr 
0 


1  2  4  n  O 

k  sec  9  e 


-2k  f  sec"  9 


-tt/z 


sin  (ak  sec  0)de 
0 


-2i  t f* 

(k^  sec^  9  +  It)  t^  e  slnh  (at  cos  0)d0  Idt 

t^  +  k  ^  sec'’’  0 
0 

The  second  Integral  of  the  right  hand  side  can  be  shown  to  vanish 
because  of  the  odd  function,  slnh. 

Hence 


00  TT 

+  I  I 

0  -IT 


2v 


=  I 
0 


v/s 


-2  k  sec" 
o 


sec^  0  sin  (ak  sec  0)  d0 
'  0 


=  R1 


[ak  ,k  ,4 
L  0  0  J 


[12] 
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If  we  expand  sin  (ak^j  sec  9)  into  a  series  of  ak^  sec  9, 


Rl[ak  ,k  j4]  =  /  e 
°  °  0 


v/2  -2k  sec^  9  n  (ak  ) 


2n  -f  1 


(-1) 


(2n  +  1) 


n=0 


X  sec^*^  ^  ^  9  d  9 


[13] 


To  evaluate  the  Integral 


tt/s  -2k  sec^  9  2v  +  1 
I  =  /  e  °  sec  9  d  9 

^  0 


[14] 


u 

we  put  sec  9  =  cosh  and  obtain 


-k 


e 

2 


00  -k  cosh  u 
r  o 

/  e 

0 


1  H-  cosh  u 


-h 


d  u  = 


^  G  (k  ) 
,v  +  1  v'  0 


00  -k  cosh  u  , 

G  (k  )  =  /  e  °  (1  +  cosh  u)  du  =  K  (k  )  -  K  (k  ) 

1  0  Q  o'  o  o'  0 


[15] 


where  K  ' (k  )  = 


d  K  (k J 


0  0 


0  '  0'  d  k 


=  -  K  (k  ) .  In  addition  we  can 
1  0 


prove  easily 

00  -k  cosh  u  n 

G(k)=/e°  (1+  cosh  u)  du 

n  o  Q 


=  G  ^  (k  )  -  G'  (k  ) 
n-1  0  n-1  o 


[16] 


J 


HYDRONAUTICSj  Incorporated 


-11- 


Hence  we  get 
-k 


I  =  ^ 

V  2^  ^ 


K  (k  ) 
'  o'  o 


K'  (k  ) 
0  0 


(n) 


where 


r 

K  (k  ) 
0  0 


K'  (k  ) 
0  0 


n\  K  (k  ) 

+  21  - ^ 


(n) 


mi  dK  (k  ) 

K  (k  )  -111  ■ 

0  '  0  dh 


d\  (k  ) 
n  0  0 


+  ( -l) 


dh" 


dh' 


n 


[17] 


We  now  perform  the  Integration  of  [13]  term  by  term  and  obtain 

(ak  "• 

^  (2n^  1)1  ^n  +  2 

n=0 


-k 


,  T  m,  ,  2n+l 
(-1)  (ak^) 


K  (k J-K:(k J 


n=0 


^n+2  r  o'  o^  o'  o', 
(2n+l)!  2  J(n+2) 


[18] 


This  series  is  proved  to  be  convergent  for  any  ak  and  k  >0 

o  0 

(see  appendix) .  When  ak^  is  sufficiently  small  this  series 
converges  rapidly.  In  computations,  the  following  recurrence 
relations  are  of  great  value. 


G  (h)  =  2  + 

n  h 


Vl  (>’)  -  °n-2 


[19] 


I 
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for  all  Integers  n.  This  can  be  easily  proved  by  integrating 
Equation  [15]  by  parts.  In  addition 


G  (h)  =  K  (h) 
o'  0 

G  (h)  =  K  (h)  +  K  (h)  [20] 

1  0  1 


On  inserting  the  value  of  R  obtained  by  the  foregoing  pro- 

o 

cedure  into  Equation  [10]  we  finally  obtain  the  non  dimen¬ 
sional  j  horizontal  component  of  the  force  at  the  doublet. 


6.  FORCE  AT  SOURCE  POINT 

Similarly,  the  force  at  point  (0,0, -f)  is  obtained 
from  [2]  in  dimensional  form 


0,0, -f 


k  ^  TT  0 

+  Re  /  / 

’’  -ir  0 


Ik  sec  9  e 


k(lu)  +  z-f) 


dkd0 


k-k  sec  0  -  Iv  sec  0 
o 


0,0-f 


1 

(a^+4f2)3/2 


3a^ 

(a^+4f^) 


+  Re 


k  LL  TT 
0 


-IT 


00 


/ 

0 


k^  "-"^)dkd0 

k-k  sec'^  0  -  Iv  sec  0 
o 


TT 
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ld<t> 

20 

dx  /0,0,-f 


k  m  TT  0 
Re-^/  / 


TT 


-TT  0 


ik  sec  9  e  dkd0 

k-k  sec^  9  -Iv  sec  9 
0 


m 


r 


K  (k  f)  +  K  (k  f)| 
o'  O  1  o  ' 


R  (0,0, -f) 


4  TT  p  m 


+ 


0,0, -f 


The  Integrals  are  evaluated  exactly  in  the  same  way  as  In 
Sections  3-5- 


7.  TOTAL  FORCE  AND  OPTIMUM  PARAMETERS 

The  total  force,  or  the  total  wave  resistance  Is,  in 
nondlmenslonal  form. 


R  =  R  (-a, 0,-1)  +  R  (0,0, -1) 


2b' 


K  (k J  + 


r^p;  ‘  °  ° 


1  + 


K  (k  ) 
1  o 


K  (k  )  +  K  (k J 


2F 


4  o '  O  1  '  o 


8b^ 

FfS 


Rlfak  ,k  ,4 
0  0  - 


where  Rl(ak  ,k  ,4)  is  given  by  Equation  [l8].  We  may  write 
0  o 


R  =  R 


-I-  R  +  R.  i. 
b  Int 


[21 


alternatively 
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where  the  nondlmenslonal  wave  resistance  (or  wave  resistance 
coefficient)  due  to  a  source 


R  =  -  <K  (k  )  +  K  (k  )| 


and  due  to  a  doublet 


and  the  interference  term 


R,  Ri  [ak  ,k  ,  4] 

int  6  L  o  o  J 

f 

Hence  R^  +  R_,  ,  is  the  total  Influence  of  the  bulb.  To  obtain 
b  int 

the  optimum  value  of  the  distance  a  which  makes  R  minimum  we 
notice  that  Rl(ak^jk^,4)  is  the  only  function  of  a  appearing  in 

the  right  hand  side  of  the  above  Expression  [21].  The  optimum 
value  of  a  is  readily  obtained  from  Equation  [l8]  with  the  aid 
of  a  numberlcal  -  graphical  procedure. 

The  optimum  value  of  b  was  obtained  by  differentiating  R 
with  respect  to  b,  putting  it  equal  to  zero,  and  solving  for  b. 


2f/  R1  fak  ,k  ,4 

f  I  o  o  J _ 

K  (k  )  +  1+7^  K(k  ) 

0  0  2k  ,  o 

o 


3, 


[22] 
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2/0 

The  optimum  values  of  a  and  br  ^  thus  obtained  are  plot¬ 
ted  against  in  Figure  4. 

8.  REMARKS 

With  these  optimum  values  of  a  and  b^  R/r~  is  plotted 
In  Figure  3  versus  F“ .  Here  we  can  see  that  a  considerable 
reduction  In  the  wave  resistance  Is  obtained  by  using  an  opti¬ 
mum  bulb.  The  force  at  the  point  of  the  doublet  In  the  pre¬ 
sence  of  the  source  with  optimum  values  of  a  and  b  Is  also 
plotted  In  Figure  3-  We  can  see  there  that  the  favorable  effect 
of  the  bulb  Is  due  to  the  pulling  force  (negative)  at  the 
doublet.  The  force  at  the  point  source  In  the  presence  of  the 
doublet  Is  not  shown  there,  out  obviously  it  Is  a  positive 
force,  which  Is  the  sum  of  the  absolute  value  of  force  at  the 
doublet  point  and  the  net  force  which  are  both  shown  in  Figure 
3.  This  reveals  that  the  reduction  of  force  on  the  half  body  - 
bulb  combination  is  entirely  achieved  by  the  creation  of  a 
large  bulb  thrust  and  not  by  the  reduction  of  resistance  on 
the  half-body  Itself.  In  fact,  the  wave  resistance  of  the 
half-body  by  Itself  would  appear  to  increase  In  the  presence 
of  the  bulb . 

The  wave  resistance  of  a  source  has  a  singularity  at 

f =0 .  The  wave  resistance  coefficient  C  Increases  as  f/r  de- 

w  ' 

creases  by  the  factor  (f/r)^  at  a  fixed  Froude  number  as  shown 
in  Figure  3  or  Equation  [10].  Although  the  shape  of  the  body 
formed  by  a  source  near  the  surface  looks  more  like  a  ship 
shape,  this  singularity  prevents  us  from  approaching  it  In  this 


way . 
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Proin  the  radius  of  the  bulb  b  shown  In  Figure  4  we 
obtain  the  strength  of  the  doublet  p  =  V  b^/2 .  This  relation 
is  used  only  for  approximating  the  si2.e  of  bulb.  The  doublet 
here  by  no  means  results  In  a  sphere.  There  is  not  only  the 
Influence  of  the  free  surface  but  also  that  of  the  source  to 
alter  its  shape  from  that  of  the  sphere.  The  form  of  the 
body  will  be  discussed,  in  Section  l4  later.  Although  the  dis¬ 
tance,  a,  between  the  doublet  and  the  source  seems  to  be  large 
compared  to  b  it  is  highly  probable  that  the  resulting  bulb  and 
bow  are  connected  by  faired  lines  (See  Figure  9A) . 


CASE  II 


9.  WAVE  HEIGHTS 


We  now  consider  the  same  problem  as  in  Case  I  except  that 
the  point  source  is  replaced  by  a  vertical  line  source  extend¬ 
ing  from  the  free  surface  to  depth  d.  In  order  to  find  the 
wave  resistance  it  is  now  more  convenient  to  use  Havelock's 
formula,  (Havelock,  1934)  rather  than  Logally's  theorem.  With 
the  former  method  it  is  necessary  to  determine  the  wave  height 
at  a  large  distance  down-stream  of  the  body. 

As  shown  in  Equation  [8]  the  x  component  of  perturbed 

velocity  of  the  flow  due  to  a  point  source  at  depth  f  with 

1 


strength  m  is  for  large  x 


x,y,0 


mk 


TT 


TT 


Re  / 


-TT 


00 


I 

0 


k(lco^  -  f^) 

1  k  sec  9  e _ dkd0 

k-k  sec^  0  -  Iv  sec  0 
o 


[23] 
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where  oo  =xcos0  +  y  sin  0 
1 

Hence,  If  we  consider  a  vertical  source  line  with  uniform 

strength  m  per  unit  length,  (or  m=m  d),  then  on  the  free 
1  1 

surface  we  have  for  large  x 


d<l) 

dx  x,y,0 


m  k 

1  o 

TT 


d  TT 


Re  /  / 

0  -TT 


k(laj^-fi) 

Ik  sec  0  e  dkd0  df, 

k-k  sec"  0  -  Iv  sec  0 
o 


,  -kd  klcu^ 

_ o  pe  j-  j-  (e  -1)  1  sec  0  e  dkd0 

^  -TT  0  k-k  sec*^  0  -  Iv  sec  0 

0 


Tr/a  -k  d  sec"  0 

=8mk/  (1-e  )  sec  0  cos  (k  x  sec  0) 

1  o  0  o 

X  cos  (k^  y  sin  0  sec"  0)  d0  [24] 

using  contour  Integration  as  was  done  In  [ll]  and  [12].  Hence 

the  wave  height  far  from  the  source  line  Is,  from  [5]  with 

V  =  0 
1 

m  Tr/a  -k  d  sec"0 

C=8-p/  (l-e°  )  sec  0  cos  (k  x  sec  0) 

'=>3  V  ^  '  o 

X  cos  (k  y  sin  0  sec^  0)  d0 
'  0 


[25] 


HYDRONAUTICS ,  Incorporated 


-18- 


Slmllarlyj  the  wave  height  due  to  a  doublet  at  (-ajO^-f)  is, 
for  large  x 

P,  2  tt/  -k  fsec-0 

ok  [I  V2  0 

C  = - ^ -  /  sec^  0  e  [sin(k  x  sec0)  cos  (k  a  sec0) 


+  cos  (k  X  sec0)  sin  (k  a  sec0)]  cos(k  y  sin  0  sec^0)d0  [26] 
'  0  o  '  0 

10.  HAVELOCK'S  FORMULA  AND  WAVE  RESISTANCE 

Suppose  the  wave  height  at  large  x, 
tt/z 

^  =  /  (P  sin  A  cos  B  +  P  cos  A  sin  B  +  P  cos  A  cos  B 

0  1  2  3 

+  P  sin  A  sin  B)d0  [2?] 

4 


where  A  =  k  x  sec  0.  B  =  k  y  sin  0  sec^0.  P^  is  a  function 
0  0  1 

of  0  and  other  physical  parameters.  Then  Havelock's  formula 

(1934)  for  wave  resistance  is  given  by 


1  Tr/2 

R=-j-p7rV^/  (P^  +  P^  +  P-  +  P-)  cos^  0  d  0  [28] 

^  n  1  2  3  4 


Hence  in  our  case  C  =  C  +  L  and 

s  ^b 
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tt/s 

R  =  l6pTr/ 

0 


-2k  f  sec  0  -k  d  sec  0 


sec  0  +  m  (1-e 

•  1 


COS0 


-k  f  sec^0  -k  d  sec^0 

-2k^“  m  [1  e  °  (1-e  °  )sec^  0  sin(k^  a  sec0)]d0 


=R,  +  R  /)  j 

b  is  il 


[29] 


where  the  resistance  due  to  the  bulb  is 


R  =  16  p  TT  k  /  e 


ttA  -k  f  sec^0 

>  r\ 


sec  9  dd 


-k  f 

=  4Trpe  °  k^p.^^K(kf)  + 

^  o  1  o^  o 


1  + 


2k  f 
0 


KJV)^ 


from  Equations  [l4]  -  [17]  where  K  ,  K  j  modified  Bessel  func- 

o  1 

tlons  of  the  second  kind.  The  resistance  due  to  the  line  source 
Is 

ir/a  -k  d  sec^  0 

R  =l6p7rm‘^/  (1-e  °  )“  cos  0  d0 

iis  1  Q 


=  16  p  TT  m“ 


■kodA 

1  -  kod  e  >|Kjk_^d/2)  -  Kjk^d/2); 


-k  d 

+  k  e  °  1k  (k  d)  -  K  (k  d))> 
^  0  I  1  o  o'  o  ' 


and  the  resistance  due  to  the  Interaction  between  the  bulb  and  the 
line  source  Is 
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-k  f  sec^e  -k  d  sec^ei 
=  -  32  pTTk^"  p  e  °  (l-e  °  | 

X  sec^  8  sin  (k  a  sec  0)  d0 
'  o  ' 

which  can  be  evaluated  by  a  series  expansion  as  in  Section  6. 


11.  NQNDIMENSIONAL  FORM 
If  we  put 

_  Vd  Vb^ 

=  m  —  p  = 

where  m  is  a  nondiraensional  coefficient,  we  nondimenslonallze 
with  respect  to  d,  e.g. 


C  =  R  =  - - 

”  f  P  V^d 


•2^2  d  -ygr  0 


k  d,  f 
o 


-  b  _  a 
b  =  —  ,  a  =  — 
d  d 


and  dropping  bars  for  convenience  we  obtain 


,6  -k  f 


K  (k  f)  +  1  +  K,  (k  f ) 

o'  o  2k  f  1  '  o  ' 

I  o  I  j 


2  1-k  e 


/k  i 

0 

-  "of) 

\2  1 

X  <^K  (k  )  -  K  (k  ) 

1  1  '  o'  o  '  o' 


_b!_  .tt/s  -^o^ 


R^l  =  -  8m  —  /  e 

d 


-k  sec  0 

l-e  o  sec^  0 


X  sin  (k^a  sec  0)  d0 


[35] 
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sistance  becomes 
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R,  +  R  . . 
b 


=-8m‘ 


rTr/s  -^0^ 

J  ® 

LO 


-kod  sec^sl  o  /  ^  ' 

ll-e  I  sec^S  sln(k^a  sec0)  d0 


-kof 


Kjk.f)  +  1  +  —  K,,  (k  f) 


o '  o 


2k  f 
0 


[40] 


The  graphs  of  the  optimum  value  of  a  and  the  optimum  value  of 


m  for  the  obtained  optimum  value  of  a  are  plotted  with 


V 

V 

respect  to  for  several  values  of  fj  the  depth  of  the  doublet^ 
in  Figures  6  and  7.  R/m^,  (R^  +  R^^)/m^j  and  R^^/m^  are  plotted 


is'^ 


versus  F^j  taking  the  optimum  values  of  a  and  b  for  each  P,^  in 
u  d 

Figure  5. 


13.  SHIP  SHAPES 

From  the  boundary  condition  [6],  we  can  see  that  the  free 
surface  behaves  like  a  solid  wall  when  the  velocity  V  is  small 

(if  P  <  .3j  7“^  <  .09  in  Equation  [6]).  Inui  (1957)  found  by 
Jj  K  L/ 

O 

comparing  his  experimental  results  with  theoretical  analyses 
that  the  effect  of  assuming  that  the  free  surface  acts  like  a 
solid  wallj  when  computing  the  stream  surface  generated  by  a 
given  singularity  distribution,  led  to  negligible  error  up  to  a 
Froude  number  of  F  =  0,7.  Accordingly  we  have  for  the  poten- 

Lf 

tlal  due  to  the  source  line  and  its  mirror  image: 


m  df 
1 


'\/x^+y^+(z-f ) 


=  m 

2  j. 


log 


z+d  +  Vx^+y^+(z+Tp 
z-d  +  ’\/x^+y^+(z-d  )^{ 


HyDRONAUTI.es,  Incorporated 


-23- 


.Por  the  doublet 

/ 

=  -  M-(x+a) 

\ 


( x+a ) ^+y^+ ( z+i )  ‘ 


3/2 


+ 


1 


(x+a)^+y^+(z-f )‘ 


Considering  a  sink  line  from  point  (L, 0, O)  to  point  (L, 0, -d) 
In  addition  to  the  doublet-source-line  system  we  find  for  the 
total  velocity  components. 


„  3<t> 

U  ,  .  _  +  V,  V  =  -  , 


Then  the  streamline  equations  are 


d4> 

dz 


dx  _  ^  _  dz 
u  ”  V  ”  w  ’ 


We  solve  this  by  Runge-Kutta-Gill ' s  (Runge,  l895i  Kutta,  I9OI; 
Gill,  195I;  or  see  Ralston  and  Wllf,  i960)  numerical  method  to 
obtain  the  three  dimensional  streamlines  starting  from  the  stag¬ 
nation  points.  Figures  9-11  show  the  body  streamlines  repre¬ 
sented  by  three  special  cases  of  different  singularity  distribu¬ 
tions  respectively:  vi^.  a  point  doublet,  a  source  line,  and  a 
sink  line  with  the  three  different  sets  of  parameters.  Figure 
17  shows  the  case  of  the  doublet  line  discussed  In  Case  III. 

Each  of  these  figures  shows  the  projections  of  the  several  body 
streamlines  on  the  x  y  plane  and  the  x  z  plane.  The  calcula¬ 
tions  were  performed  on  the  IBM  1620.  In  making  these  computa¬ 
tions  It  was  necessary  to  exercise  extreme  care  In  the  selection 
of  Interval  increments  at  the  neighborhood  of  the  stagnation 
point  and  in  the  vicinity  of  large  curvature.  The  resulting 
body  streamlines  are  smooth  curves  with  a  hollow  place  although 
the  distance  between  the  doublet  and  the  source  appears  to  be 
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too  large  to  be  a  single  connected  body. 

14.  DISCUSSION 

The  optimum  distance  between  the  doublet  and  the  source 
llnej  and  the  optimum  size  of  the  bulb  Increase  with  the  flow 
speed  and  depth  of  the  doublet  as  shown  In  Figures  6  and  7* 
However,  the  depth  effect  of  the  doublet  on  the  wave  resistance 
decreases  with  flow  speed  as  long  as  we  take  the  optimum  bulb 
size  and  the  optimum  distance  between  the  doublet  and  the  source 
line  at  each  selected  depth.  This  Is  shown  In  Figure  8. 

The  size  of  bulb  Is  considerably  large  for  the  reasonably 
high  Froude  number.  When  F^^  >  2  the  width  of  the  bulb  becomes 

larger  than  the  beam  length  as  shown  In  Figures  7  and  10. 

Hence  the  bulb  may  be  considered  rather  as  a  blunt  bow  of  ship 
than  an  appendage.  However,  as  we  notice  In  Figures  10  -  11, 
the  hollow  place  between  the  bulb  and  ship  body  may  contribute 
to  the  flow  separation. 

Assuming  no  serious  separation  takes  place  the  reduction 
of  the  wave  resistance  due  to  the  bulb  Is  remarkably  great. 

Figure  5  shows  that  the  effect  of  an  optimum  bulb  is  to  reduce 
the  wave  resistance  by  more  than  60  percent  of  that  contributed 
by  the  line  source  alone  at  all  Froude  numbers. 

The  wave  resistance  coefficient  of  the  source  line  bow 
alone  Is  also  shown  In  Figure  5.  It  has  a  Jump  at  F^=  0,  and 
It  decreases  as  F^increases.  Since  this  is  a  nondlmenslonal 
quantity  It  does  not  mean  that  the  wave  resistance  Itself  has 
a  jump  at  =  0, 
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CASE  III 

15.  WAVE  RESISTANCE  OF  THE  SYSTEM  OF  A  DOUBLET  LINE  AND  A  SOURCE 
LINE 


As  shown  in  Case  II  the  optimum  size  of  the  bulb  is  re¬ 
markably  largej  especially  for  higher  Proude  numbers.  We  will 
now  Investigate  the  result  when  the  doublet  strength  is  dis¬ 
tributed  linearly  along  the  vertical  line  from  (-a,Oj~f^)j,  to 

(-ajOj-fg)  in  front  of  the  same  source  line  as  in  Case  II. 

If  we  take  [i  =  +  fp^  for  the  doublet  strength  per  unit 

length,  we  have  for  the  wave  height  far  behind  the  system  due 
to  the  bulb,  from  Equation  [26] 


C 


b 


8k  2 
_ 0_ 

V 


jrA 

0 


-k^f  sec^0 
e 


sec'^  e 


X  [sin  (k^x  sec0)  cos  (k^a  secS)  +  cos  (k^x  seed) 
X  sin  (k^a  secS)]  cos  (k^y  sin  9  sec^B)  d  f  d  0 


8  .f/s 
^  0 


-k  f  sec  0 

k^sec^0^e  °  ^  (p^fi+M-i  )“  o 


■kQfgSec  0 


X  (^^_f2+^il  )f+ 


-k^f-  3ec^0  k^f^sec^0 
0  1  o  a 

;  -e 


X  [sin  (k^x  sec0)  cos(k^a  sec0)  +  cos  (k^x  sec0) 


X  sin  (k^a  sec0)]  cos  (k^y  sin  0  sec^0)  d0  [4l] 


and  due  to  the  source  line,  from  Equation  [25] 
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=s  V 


ir/z  ~^o' 

'  1  -  e 


-k„d  sec  0 


sec  0  cos  [k  X  sec  0] 
o 


X  cos  (k  y  sin  0  sec^  0)  d0 
'  o 

Using  Havelock’s  formula  (27)  and  (28 )j  we  have  for  the  total 
wave  resistance^  due  to  the  doublet  and  source  lines 


ttA  ^  -l^^d  sec  0 

R  =  l6  Trp  /  ~  ®  ° 

0  ^ 


^  sec  0  ,2 

+  sec  0<e  (p  f^+p.^  )  +  e 


2  -Bk^f^sec  0 


(u  f  +\i  ) 

Q  O  J- 


-k^(f, +f_ )sec^0 

-2  e  I'*' 

-2k  fisec^0  -2k  fgSec^© 

X  e  °  )p-2+  e  -’^2 


-kQ(fj^+f2  )sec^0 


cos  0 


-2k  f  sec^0  -2k„fpSec^0 

+  p  ^  e  cos^0+p  ^  e  cos®0-2|j,  ^ 

2  2  2 


-kQ(f^+f2  )sec^0 


^0  „  [  -knfisec^© 

cos®0+  -2k  m  <e 
o  1  1 


-k  fgsec  0  -ko(fj_+d)sec  0 


(f^Pg+p^ ) 


-k-(f  +d)  -k„f,  sec^0  f  sec  0 

i  ^  (fgli^+M-^^  )>-2ra^|j,^cos®0  e  °  -e 

-k^(f. +d)sec^0  -k  (f  +d)sec®0| 

I  c  ^  +  e  sin  (k  a  sec  0)^^d0[42] 
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16.  NONDIMENSIONAL  FORM 


We  perform  the  integrations  using  the  methods  taken  in 
Equations  [11]^  [15]j  and  [l6],  and  nondlmenslonalize  physical 
variables  with  respect  to  m^^  and  d  as  in  Case  II,  i.e. 


„  -  R  —  -  ^2 

C  =  R  =  -  ,  11,  =  — ^  ,11  =  —  ,  m,  =  — r— 

w  T  2-,r2^s  ^  m  d  *^2  m,  i  4 

—  p  m  V  d”^  11  1 


=  -l=r  ,  k  =  k  d,  f  =  f ,  /d  ,  f  =  f/d 

and  dropping  bars  for  convenience  we  obtain 

R  =  R,  +  R«t^  +  R/)». 
is  ib  iii 

r  \  /k  I  /  _v  I 

R,  =  2  l-k<^K  -  K  -^le  +  k  e  °<^K  (k  )-K  (k  V 

is  [  oj  1  2  /  o  2  I  0  I  1  ^  0-^  o'  0*^1 

•w  J 


''ib  4  2 


-k„f 


2  -^0^2 

K  (k  f  )(u  f  +il,  )  +  - 

o  '  O  1  2  1  *^1  '  2 


X  K  (k  f„ )  (p,  f  +11  )  -e 
0  '  0  2  '  '■^2  2  ^1  ' 


+  k  f  e~^o^i  K,  (k  f  )-K  (k  f  )  (p,  f  +p^  )p 
„  2  0  1  1'  o  1'  o'  o  1'  '  1  1  I'^z 

P  L  ' 

d 


(p^fi+Pi)(p2f2+M.,) 


-k  f. 


+  "o^2  ^  °  (^1(^2)  -  ("2^2  ^^2 

)  -k  (f  +f„)/2| 

-  ^  e  °  ^  ^  <K  k  (f  +f  )/2  -K  k  (f  +f  )/2 

2  )  1  o '  1  2 '''  0  0 '  1  2 


4  of 

^  ,.2  „  o  1 


2p,P,+P3^(f^+fJ  +  3  P^  e  o  ^  2k^f^K^(k^fJ-(2k^f  -l)K^(k^fJ  k^f^ 


k  f  <2k  f  K  (k  f  )-(2k^f  -1)K  (k  f  ) 

02]  020'02'  '  Og  ^I'Os^ 


X  K  (k  (f  +f  )/2  -Ik  (f  +f  )-  1)  K  k  (f  +f  )/2 

o|o  i  2'^  \o  i  2  /  i|o'i  2'' 
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/  r  ^  ^  1  r 

^  ,  0  -  R1  ak^, 

\  L  J  L 


ko(Vl) 


r  k  f  -]  r  k  (f^+l) 

,  0  -  Ri  ak^,  -2-2 -  .0 


r  r  ^  ^  1  r  ^  ^ 

4  J  R1  ak^,  -  2  -  R1  ak^,  -|^  ,  -  2 


r  k  (f  +1)  -)  r  k  (f  +1) 

R1  ak^,  -2~ -  ,  -  2  +  Ri  ak^,  ,  -  2 


where  Rl  [gjk,  n]  =5  / 


-2k  sec  6  n 


sec  9  sin  (g  sec  0)  d0 


which  Is  defined  In  Expression  [12]j  and  Is  evaluated  In  Equa¬ 
tions  [12]  -  [18]  by  the  series  expansion  In  g. 

1'7.  OPTIMUM  PARAMETERS 

The  optimum  values  of  u  and  u  which  make  R.,  -i-  R... 

^2  -gb  ££,1 

minimum  may  be  obtained  by  the  usual  methods.  We  differentiate 
R  partially  with  respect  to  and  p^  respectively^  put  the  re¬ 
sults  equal  to  zeroj  and  solve  the  resulting  two  simultaneous 
equations  for  the  optimum  values  of  p^  and  p^.  The  equations 

for  p^  and  p^  resulting  from  the  above  mentioned  differentia¬ 
tion  are 


U  = 
1 


Y  -  p  X 
1  2 


Y3-  P,X 


IJ-  = 
2 
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3^d 


2  -  O  1 


f  e 


K  (k  f  )  +  f  ^  e 

o '  0  1  '  p 


k  f 

2  °  ^  K  (k  f  ) 

P  0^02' 


+ 


K 


3P. 


+ 


3P. 


/f " 

-k  f 

+  f 
^  1 

e  °  "  + 

i 

2 

X  e  K,  (k  )  - 

1  '  o  1 


(f  +f  ) 

0  1'  1  --  O  '  1  P 


X  K, 


0  '  1  2  ' 


2P. 


■1 1  2 


The  optimum  value  of  the  distance,  a  which  makes  minimum 

is  obtained  numerically  as  in  Cases  I  and  II.  The  optimum  values 

of  u  and  u,  were  calculated  for  various  values  of  f  and  f  . 
12  1  2 

The  results  of  these  calculations  are  shown  In  Figures  12  and 

13.  The  wave  resistance  coefficient  for  the  obtained  optimum 

parameters  for  each  Proude  number  is  also  calculated,  but  this 

Is  the  same  as  that  of  Case  II. 


18.  EFFECT  OF  STERN 

The  stern  can  be  represented  as  a  uniform  sink  line  from 
(L, 0, 0)  to  (L, 0, -l)  whose  total  strength  Is  the  same  as  the 
source  line.  Here  L  Is  nondimenslonalized  with  respect  to  d. 

The  wave  resistance  due  to  the  source  and  sink  lines  alone  ex¬ 
cluding  Interferences  Is  as  twice  as  that  due  to  the  source  line 
alone.  However,  the  effect  of  the  Interference  between  the 


HYDRONAUTICS,  Incorporated 


-31- 

singularlties  representing  the  forward  and  aft  part  of  the  ship 

may  be  extrdmely  important.  The  wave  resistance  between  the 

source  and  sink  lines  R  . ,  is  evaluated  using  the  method  of  sta¬ 
ll 

tionary  phase  (see  e.g.  Stoker^  1957)  In  nondlmensional  form. 


R  .  =  -  4  / 
0 


Tr/2 


1  -  e 


-k  sec  e 


cos  9  cos  (k  L  sec  0)  d0 
'  o  ' 


2ir 

|2 

(  -^0 

k  L 
o 

i'  ■ 

cos  (k^L  +  f) 


[48] 


which  is  approximately  valid  when 


<  .4  (Inui,  1955)  ■ 


Similarly  the  wave  resistance  due  to  the  interference 
between  the  doublet  line  and  the  sink  line  is 


+  H, 


1-  e 


-k^ 


p-T'-je  ° 


e-‘'oA 


rik)\  sinjkJa+L)+f 


[‘>9] 


The  optimum  values  of  and  \i^  change  very  slightly  when  the 

sink  line  stern  is  considered  (see  Table  l). 

The  total  resistance  for  the  system  including  the  sink 
line  stern  and  optimum  doublet  lines  is  plotted  in  Figure  15  for 
Proude  numbers,  F  =  0.2,  0.25,  and  0.3.  The  numerical  values 

ij 

of  the  stern  effect  for  each  Proude  number  are  also  shown  in 
Table  1.  Here  we  can  see  that  the  interference  between  the 
stern  and  the  bow  plus  bulb  is  in  general  much  less  than  that 
of  the  stern  and  the  bow  alone.  That  is,  the  bulb  has  an  effect 
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of  smoothing  out  the  humps  and  hollows  of  the  resistance  curve 
to  a  considerable  extent.  Besides j  the  magnitude  of  the  inter¬ 
ference  resistance  between  the  bulb  and  the  bow  is  much  larger 
than  either  the  Interference  resistance  between  the  stern  and 
the  bow  or  between  the  stern  and  the  bulb. 

19.  DISCUSSION 

Under  the  optimum  conditions^  the  effect  of  the  line 
doublet  is  similar  to  that  of  the  point  doublet  of  equal  total 
strength.  The  optimum  position  of  the  center  of  gravity  of  the 
line  doublet  Increases  in  depth  and  mo  es  forward  with  increas¬ 
ing  Froude  number  while  at  the  same  time  Increasing  in  total 
strength.  This  is  shown  by  Table  1,  Figures  12  and  13.  The 
total  effect  of  the  line  doublet  on  the  wave  resistance  is  al¬ 
most  exactly  the  same  as  that  of  the  point  doublet. 

The  wave  resistance  curves  due  to  the  systems  of  source 
line  and  the  doublet  line  alone  optimized  for  Froude  numbers 
Fj^  =  0.2j  0.25j  and  0.3  are  shown  in  Figure  lA.  The  wave  re¬ 
sistance  at  a  little  lower  Froude  numbers  than  the  optimum 
Froude  numbers  is  larger  than  that  due  to  the  source  line  alone, 
but  at  Froude  numbers  in  the  vicinity  of  and  greater  than  that 
for  which  the  line  doublet  was  optimized,  the  effect  of  the  line 
doublet  (bulb)  is  always  to  reduce  the  wave  resistance. 

The  variation  of  the  strength  of  the  doublet  line  (assumed 
to  be  linear)  is  very  small  in  general  unless  we  make  the 
doublet  line  long  enough  so  that  its  upper  end  is  very  near  to 
the  free  surface.  In  this  case  the  optimum  slope  is  quite  large 
for  low  Froude  numbers .  However,  the  upper  end  of  the  linear 
doublet  line  should  not  be  too  close  to  the  free  surface  for  op¬ 
timum  interference.  In  fact,  it  can  be  seen  from  Table  1  that 
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the  optimum  depth  of  the  upper  end  increases  with  Increasing 
Proud e  number. 


2  0.  APPLICATION 

If  we  want  to  Improve  the  ship  shape  of  length  L,  whose 
waterline  Is  given  by  the  polynomial  (G.  P,  Welnblum,  1950)  as 

A 


n 

n+l'  ^ 


n=0 


then  the  singularity  distribution  in  the  sense  of  Michell's 
ship  is 

5 


n 


A  X 
n 


n=0 


The  height  of  bow  wave  at  a  large  x  is  given  by  (see 
Equations  [25]^  [26]), 


r  = 

^bs  Vk  ■' 


o  0 


1-  e 


-k^z^sec  0 


21  a 


4i  a 


+ 


4 


°  k  ^sec^0  k  '^sec‘^0 
0  0 


X  sin  (k  X  sec  0)  + 


3J  a. 


5i  a. 


k  sec0 
o 


k  ®sec®0  k  ^sec^0 


cos(k^x  sec9))> 


X  cos  (k^y  sin  0  sec^  0)  d0 


[50] 


If  we  combine  the  bulb  whose  wave  height  at  a  large  x  is 


tt/  *4:„f  sec^0 

C.  =  -  -f-  sec^0  e  ° 


<sin  (k  X  sec0)  cos  (k  a  sec0) 


0 


+  cos  (k^x  sec0)  sin  (k^a  seed)  cos  (k^y  sin0  sec^0)^  d0 


the  wave  resistance  due  to  this  bulb  and  bow  are  obtained  from 
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Equatlon  [29]  and  are  given  by 


R 


=  l6p7r  '  !l- 

0  ' 


-k  z  sec^0'  I  ® 


O  ^2 


4!  a 


4 


+ 


^0  k  ^sec^S  k  ^sec^Ol 
0  0 


-k  f  sec^d 

-  pk  ^  e  sec  0  cos  (k  a  sec  0) 


cos^  0 


+  cos^0 


1- 

2 

f  sec  9\ 

(  n  3-  \  5!  a. 

r  J 

L  -j  - 

. 

i  k  ^sec0  k  ^sec^0  k  ^sec^0  ■ 

-  pk  ^  e  ° 
0 


-k  f  sec  0  T 

sec^0  sin  (k^a  sec  0)j 


d0 


All  .the  Integrals  Invol/ed  here  are  of  the  form 

I  =  g  ®  ggQ  0(3g  (shown  in  Equation  [14]) 

0 

if  we  consider  the  expansion  of  cos  (k^a  sec0)  and  sin  (k^a  sec0). 

Hence  the  wave  resistance  in  this  case  may  be  evaluated  exactly 
in  the  same  manner  as  in  Section  5- 

The  optimum  parameters,  f,  p,  and  a  may  be  determined  by 
the  methods  used  In  Sections  7  and  12. 

If  we  Include  the  stern  the  expression  for  the  wave  re¬ 
sistance  becomes  a  little  complicated.  However,  the  forms  of 
the  Integrals  are  the  same  as  those  encountered  previously  when 
the  method  of  stationary  phase  was  used  for  F  <  0.4. 

Since  the  effect  of  the  bulb  Is  largely  due  to  the  inter¬ 
ference  of  the  bulb  with  the  bow  wave  (Takahel,  i960),  we  need 
not  be  concerned  about  the  Influence  of  the  stern  In  dealing 
with  the  determination  of  the  optimum  bulb  parameters  even  though 
they  vary  slightly  as  shown  in  Section  I8. 
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The  elementary  bow  waves  (Havelock,  193^a)  of  usual  ships 
are  combinations  of  sine  waves  and  cosine  waves  as  shown  in 
Equation  [50].  Since  the  best  position  of  the  doublet  for  the 
sine  bow  waves  can  be  easily  shown  to  be  at  the  bow  itself,  the 
distance  between  the  doublet  and  the  bow  in  the  case  of  usual 
ships  will  not  be  as  large  as  is  in  the  case  of  a  source  line 
ship.  Hence  the  series  expansion  used  in  finding  the  distance 
may  be  evaluated  more  easily. 
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APPENDIX 


Consider 


2 

TT  /  n  “^Nsec0  pv  .  » 

R1  (tjhj2v)  =  /  e  sec  Q  sin  (t  sec  0)  d0. 


Expanding  sin  (t  sec  6)  we  have  for  the  integrand 


-l)  (t  sec  0)  ®  2V  ^  „ 

' - - r-v — ,  -  e  sec  0  +  R  ^ 

(2  n+l)  m 


,  ,2m+3  ,  2- 

.  I  (t  sec0)  -2h  sec  0  av 

where  R^  is  the  remainder  and  <  - ( 2m+3 )'  !  ^ 


sec^  0  [Al] 


in  0  <  0  <  Tr/2  by  Taylor's  remainder  theorem.  It  then  follows 


that 


2  n+l  I 


Ri  (tjh,2V)  =  - 


n-n  (2n+l);2 


-  K  (h)-K'  (h)'-  +  /  R  d0 

n+v  '  0^  ■  A  ra 

(n+v ) 


To  prove 


Ri  (t,h,2v)  = 


- ^ -  K  (h)-K'(h) 

L  J(n+v) 


we  have  only  to  prove 


R  d0  =  0  . 
m  -*■  00  j  m 


By  the  inequality  [Al], 

/  /  -Ph  Qor.2o  ,2m+3  2V+2m+3. 

Va  r,  ^  t  sec  0 

j  R  d0<  , - 

J  m  -  i  (2m+3)  ; 


K  (h)-K'(h)\ 

1  o'  '  o'  '  , (m+i  +  v, 


(2m+3).'  2 


m+i  +  v 
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Nextj  for  later  use  we  have  to  prove  an  inequality 


v-i 

K  (h)  <  ^ -  vl  K  (h) 

V  —  .  V-l  1 

h 

where  0  <  h  <  1  and  v  is  a  positive  integer. 
From  the  recurrence  formulae 


[A2] 


Hence  by  the  mathematical  induction,  the  inequality  [A2]  is  true. 
From  the  recurrence  formula 

-  K  '  (h)  =  ^  K  (h)  +  K  (h) 

V  '  '  2  \  v-l  '  '  V+l  '  'I 

we  can  prove 
V  (h) 

(-1)  -  <  K  (h)  +  K  _  (h)  +  ...  +  K  (h)  when  v  is  odd 

dh^  ^ 

<  K^(h)  +  2^^^  +  ...  +  when  v  is  even 
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By  definition 


K  (h)-K'(h).  . 
0^  '  0^  '(n) 


K  (h)-(  ) 

o'  '  'I''  dh 


dK  (h)  ^  d^K  fh) 


Q 


dh‘ 


n  dV  (h) 


+  (-1) 


dh 


n 


[A5] 


Using  the  Inequalities  [A2]  and  [A4]  In  Equation  [A5]  we  can  prove 


n-i 


^Ko(h) 


Kj(h)|>  <  -  K^(h)  for  V  <  h  <  1 

{n) 


R  d0 


< 


.am+3  m-raJKj^(h) 


h"’  2®''‘''(2ra+3): 


The  left  hand  side  will  approach  zero  when  m  becomes  larger  than 
t^/h.  When  h  >  1  the  proof  Is  simpler. 
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source  line,  and  the  Interference  defined  by  Equations  [45], 
[^4],  and  [46]  respectively 

R  =  R.  +  R«.  +  R«/7. 

Us  ib  iii 


R, .,R  .:  Wave  resistance  Interference  between  the  source  and 
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sink  lines,  and  between  the  doublet  and  sink  lines  respec¬ 
tively 

R^:  Total  wave  resistance  Including  the  stern  with  p  and  p 
t  °  20 
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HYDRONAUTICS,  Incorporated 


TABLE  lb. 

OPTIMUM  PARAMETERS  AND  THE  WAVE  RESISTANCE  (d/L  ^  O.O7) 


''l 

.35 

.4 

.45 

.5 

^^10 

.28731 

.32629 

.36552 

1.4464 

-1.7513-05 

-2.1465-05 

-1.7576-05 

2.8336-07 

f 

1 

.04 

.04 

.04 

.06 

.07 

.07 

.07 

.07 

a 

,1145 

.127 

.1405 

.15 

^ib 

8.4268 

6.5796 

5.0410 

3.9724 

13.043 

10.311 

8.2791 

6.7546 

^iil 

-16.853 

-13.159 

-10.082 

-7.9447 

R 

4.2683 

3.7315 

3.2381 

2.7823 

’’u 

4.7276 

-2.9417 

-2.6130 

-.17448 

".I 

-1.8833 

2.9071 

.26984 

-1 . 0252 

’’t 

20.160 

14.008 

9.1740 

8.3372 

.31814 

.25420 

.35574 

1.6330 

^23 

3.8745-05 

-9.3216-05 

-2.3353-05 

7.9889-07 

p^o^M-go*  Values  of  and  in  the  negative  optimum  doublet 
strength  p  =  +  fp^j  without  considering  the  stern. 


7.5706-02  =  7. 5706- 10'^ 


Wave  resistance  due  to  the  doublet  line,  the 


source  line,  and  the  Interference  defined  by  Equations  [45], 
[44],  and  [46]  respectively 


R  =  R/i_  +  R/11-  +  R 


is 


ib 


7il 


R, ,,R  Wave  resistance  interference  between  the  source  and 
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tively 
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FIGURE  2-  CONTOUR  OF  INTEGRATION 


FIGURE  3  WAVE  RESISTANCE  OF  HALF  BODY  WITH  OPTMIZED  BULB 
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FIGURE  5- WAVE  RESISTANCE  OF  A  SOURCE  LINE  AND  OPTIMUM  DOUBLET  LINE  AT  EACH  FROUDE  NUMBER 
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FIGURE  8-  OPTIMUM  DEPTH  OF  DOUBLET 
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FIGURE  15-WAVE  RESISTANCE  FOR  SOURCE  LINE  SHIP  INCLUDING  STERN 

(Sink  Line) 
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FIGURE  16- WAVE  RESISTANCE  INTERFERENCE  WITH  STERN  (Sink  Line) 
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FIGURE  17-  BODY  STREAM-LINE  SHAPE  DUE  TO  A  DOUBLET  LINE  AND  A  SOURCE  LINE 
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